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Abstract 

Background: Zearalenone is a nnycotoxin produced by several species of Fusarium genus, nnost notably Fusarium 
graminearum and Fusarium culmorum. This resorcylic acid lactone is one of the most important toxins causing 
serious animal and human diseases. For over two decades it has been known that the mycoparasitic fungus 
Clonostachys rosea (synonym: Gliocladium roseum, teleomorph: Bionectria ochroleuca) can detoxify zearalenone, 
however no such attributes have been described within the Trichoderma genus. 

Results: We screened for the presence of zearalenone lactonohydrolase homologs in isolates of Clonostachys and 
Trichoderma genera. We report first finding of expressed zearalenone lactonohydrolase in Trichoderma aggressivum. 
For three isolates (7. aggressivum, C. rosea and Clonostachys catenulatum isolates), we were able to reconstruct full 
coding sequence and verif/ the biotransformation ability potential. Additionally, we assessed progression of the 
detoxification process (in terms of transcript accumulation and mycotoxin decomposition in vitro). 
In silico, search for origins of zearalenone lactonohydrolase activity in model fungal and bacterial genomes has 
shown that zearalenone lactonohydrolase homologs form a monophyletic fungal clade among the a/b hydrolase 
superfamily representatives. We corroborated the finding of functional enzyme homologs by investigating the 
functional sites (active site pocket with postulated, noncanonical Ser-Glu-His catalytic triad) conserved in both 
multiple sequence alignment and in homology-based structural models. 

Conclusions: Our research shows the first finding of a functional zearalenone lactonohydrolase in mycoparasitic 
Trichoderma aggressivum (an activity earlier characterised in the Clonostachys rosea strains). The supporting evidence 
for presence and activity of functional enzyme homologs is based on the chemical analyses, gene expression 
patterns, homology models showing conservation of key structural features and marked reduction of zearalenone 
content in cultured samples (containing both medium and mycelium). Our findings also show divergent strategies 
of zearalenone biotransformation ability (rapid induced expression and detoxification vs. gradual detoxification) 
present in several members of Hypocreales order [Trichoderma and Clonostachys genera). The potential for 
lactonhydrolase activity directed towards zearalenone and/or similar compounds is likely ancient, with homologs 
present in several divergent filamentous fungi among both Sordariomycetes {Bionectria sp., Trichoderma sp., 
Apiospora montagnei) and Leotiomycetes {Marssonina brunnea f. sp. 'multigermtubi'). 

Keywords: Zearalenone lactonohydrolase, Zearalenone, Clonostachys sp, Trichoderma sp, Fusarium sp, Mycotoxin s. 
Epoxide hydrolase. Homology modelling 
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Background 

The resorcylic acid lactones are mainly produced by 
fungi belonging to Hypocreales order (e.g. F, grami- 
nearum, Hypomyces subiculosus, Pochonia chlamydos- 
poria). Majority of the known compounds is bioactive 
[1]. The most widespread (due to its potential for accu- 
mulation in food and feed) is zearalenone (6-(10-hy- 
droxy-6-oxo-trans-l-undecenyil) -resorcylic acid lactone). 

Zearalenone (ZEN) - a mycotoxin produced by several 
species of Fusarium, most notably F, graminearum and 
F, culmorum - has relatively low acute toxicity, but it ex- 
hibits distinct estrogenic and anabolic properties [2], due 
to its ability to couple with the estrogen receptor. The 
examination of ZEN metabolism in swine, rat and chick- 
ens liver resulted in a- and |3-zearalenol as products (the 
latter compound being recognized as non toxic) with in- 
dication that a-zearalenol binds to estrogen receptors 
10-20 times stronger than ZEN and 100 times stronger 
than |3-zearalenol. The long-term effects of ZEN expos- 
ure include genotoxic and carcinogenic effects e.g. [3,4], 
as well as variety of reproductive disorders in animals 
e.g. [5-7]. 

In vivo, zearalenone has been proven to exhibit signifi- 
cant fungistatic effects and is thought to contribute one 
of the key mechanisms of competition between producer 
and non-producer species [8]. In keeping with this, abil- 
ity to detoxif)^ zearalenone is thought to confer a consid- 
erable adaptive advantage to competing fungal taxa [9]. 

Among the fungi of Hypocreales order, the mycopara- 
sitic fungus C. rosea was long known to degrade zearale- 
none [10]. The exact mechanism of detoxification was 
determined in form of zearalenone-specific lactonase 
(zearalenone lactonohydrolase) enzyme {zhdlOl) which 
catalyzes the hydrolysis of ZEN, a process followed by 
spontaneous decarboxylation [11]. The end products ex- 
hibit both significantly lessened toxic effects and a de- 
creased affinity for estrogen receptors. 

To this day, independent detoxification mechanisms 
have been reported both in fungi {Trichosporon mycotox- 
inivorans) [12] and in bacteria {Rhodococcus pyridinivor- 
ans) [13]. However, a systematic screening of potential 
biocontrol agents (divergent fungi of Hypocreales order - 
mainly Clonostachys sp. and Trichoderma sp.) for lacto- 
nohydrolase activity and expression patterns has not, to 
our knowledge, been described in literature. 

In this study, we present the results of screening a 
combined collection of Trichoderma and Clonostachys 
isolates, for strains with functional lactonohydrolase 
homologs and confirmed biotransformation ability. We 
report the first finding of a functional zearalenone lacto- 
nohydrolase in T, aggressivum. We also present results 
of an inquiry into the evolutionary basis of potential 
resorcyclic acid lactonohydrolase activity in filamentous 
fungi. 



Results 

Population screening for potential biocontrol agents 

Taxonomic identification of isolates used in the screening 
was carried out with use of both morphological (mycelium 
and conidia morphology) and molecular techniques (ITS 
and TEF sequences; Th2/Th4 marker [14]). 

We found seven pairs of primers amplif)^ing overlap- 
ping products nested within the zearalenone lactonohy- 
drolase coding sequence (products of ca. 300 bp). Total 
of seventy nine isolates belonging to the Trichoderma 
and Clonostachys genera were tested for the presence of 
the gene. For three isolates (C. catenulatum - AN 169, 
C. rosea - AN 154 and T. aggressivum - AN 171) we ob- 
tained full length products for all seven primer pairs, 
which allowed for the assembly of complete coding se- 
quences for respective lactonohydrolase homologues 
(see Figure 1 for multiple alignments with model Clonos- 
tachys sequence as well as reference model species se- 
quences). The three isolates were further investigated in 
detail. GenBank accession numbers: AN 169 - KF 
515222, AN 154 - KF 515223, AN 171 - KF 515221. 

Verification of biotransformation ability potential in 
isolates of Clonostachys sp. and isolate 
of Trichoderma sp. 

The fastest mycotoxin decomposition was observed in the 
isolate AN 169 (C. catenulatum), where after 24 hours the 
levels of ZEN were found to have declined below detectable 
levels (complete biotransformation ability). In the other two 
cases, the process progressed much slower. In case of iso- 
late AN 154 (C rosea), two days after incubation the 
concentration of ZEN decreased below 50% of initial 
concentration. In AN 171 culture {J. aggressivum) compar- 
able level was achieved after six additional days. In both 
cases, after full eight days of incubation the concentration 
of ZEN in the medium dropped by approximately 80-90% 
(see Figure 2). 

Zearalenone lactonohydrolase gene expression in isolates 
of Clonostachys sp. and isolate of Trichoderma sp. 

Expression of zearalenone lactonohydrolase gene was tested 
via quantitative RT-PCR (with |3-tubulin as reference gene). 
The isolate AN 171 (J, aggressivum) isolate exhibited over 
16-fold induced increase in zhdIOl expression 2 hours after 
zearalenone exposure (which corresponds with results of 
chemical analysis showing gradually expressed biotrans- 
formation ability potential). Conversely, the two other 
isolates AN 154 (C rosea) and AN 169 (C. catenulatum) 
exhibited different expression patterns. The AN 169 isolate 
(the most effective detoxifier) accumulates higher transcript 
levels slowly but consistently over the period of days, while 
AN 154 most likely presents constitutive varying en- 
zyme activity (as evidenced by low slope/plateaus in 



Popiel et al. BMC Microbiology 201 4, 14:82 Page 3 of 1 2 

http://www.bionnedcentral.conn/1471 -21 80/1 4/82 



gi|2ie93123|dbj|AB076037.1| | 

AN171 L 

AN169 ||f|G||| 

AN154 r 



|G||||G|G| II 
■GlIllGlGl 11 
IGIIIIGIGI - 



GGlllTGGGi |G|GG|GT|G 

GGlllTGGGi |G|GG|GT|G 

_ GGlllTGGGi iGlGGlGTiG 

AN154 GT|||T|T|| GG|||TGGG| |G|GG|GT|G 




TTllllGlll 




limiffGllGl lllllliGllG 97 

iimnGiiGi iiiinRGiiG ioo 

TTllllGlll GlIlllGlll limnGllGl |||||inG|iG 100 
GlIlllGlll lllf f GllGl llllf f Gif G 100 



GlIlTGGlGl 
GlIlTGGlGl 
G|||TGG|G| 



TiGllGGGll 
TTGllGGGll 
TTGllGGGll 
TTGllGGGll 



TTTlTTTGTG G|GGTTTT|T | 
TTTlTTTGTG G|GGTTTT|T L 
TTTlTTTGTG G|GGTTTT|T ITgI 
TTTlTTTGTG G|GGTTTT|T T ' 



G||TG|||GT GGIITTGIIG 197 

G||TG|||GT GG||TTG||G 200 

GllTGlllGT GGllTTGllG 200 

G||TG|||GT GGllTTGllG 200 



gi|216g3123|dbj|AB076037.1| GG|TGG||G| ||T||T||G| |G||GTG||T |GTGTTT||| 

AN171 GGlTGGllGl ||T||T||G| |G||GTG||T |GTGTTT||| 

AN169 GG|TGG||G| ||T||T||G| |G||GTG||T |GTGTTT||| 

AN154 GGlTGGllGl ||T||T||G| |G||GTG||T |GTGTTT||| 



IIGIIGGIII 
llGllGGlll 



IGIIGGI 



GGlGGTTlll |G|||TT|GG |||T|||TTT T||TG||GG| TT|G||TGG| 297 

GGlGGTTlll IgIIITTIgG |||T|||TTT T||TG||GG| TTIgIITGGT 300 

GGlGGTTlll |G|||TT|GG |||T|||TTT T||TG||GG| TT|G||TGGT 300 

GGlGGTTlll iGlllTTiGG |||T|||TTT T||TG||GTG |T||TT|GGT 300 



gi|21693123|db)|AB076037.1| |G|6||G|G!I |||TG|T|G| -TT|T||||| G| -|G|||GG |GG||T|TT - |G|TTGT|T| TGT|G||T|T G|T|T||||| TTG|GG|||| T|||||T||| 394 

AN171 IGIGIIGIGK G||TG|T|G| |TT|T||||| GT||G|||GG IGGIITITII I|TTTGT|T| TGT|G||T|T GiTlTlllll TTG|GG|||| TlllGlTllT 400 

AN169 IGIGIIGIGK G||TG|T|G| |TT|T||||| GT||G|||GG iGGlltilftll ||TTTGT|T| TGT|G||T|T G|T|T|mi TTG|GG|||| P 

AN154 |||G||G|G|I |||TG|T|G| -TT|T||||| G| |G|||GG |GG --|1|| IITTTGTITI TGT|G||T|T G|T|T|mi TTG|GG|||| | 



gi|2ie93123Idbi|AB076037.1| | 

AN171 |TT|G|G||| 
|TTiG|G||| 

AN154 r 




llillll IIIIIIIGIT ITIIITTGTI TTGi|||||| | 

::i::::git itiiittgti ttgiiiiiii i 

iiiiiiigit itiiittgti ttg||||||| | 

IIIIIIIGII ITIIITTGTI TTGVIflllll | 



I llGllllGll IGIIIIIITG iGlllTlGll G|||||T|T| 494 

ITG IGIITTIGII GIIIIITITI 500 

ITG IGIITTIGII GIIIIITITI 500 

IIGIIIIGII IGIIIIIITG IGIITTIGII G|||||T|T| 495 



GGllllllGl GllGGlGGGT ||TGG||||g ||gttgt||t |gtt|||g|t gg||t|gg|g IItgIIIgIt 

GGllilllin GGnilffillGl GllGGlGGGT |||GG||||G ||GTTGT||T |GT|||TG|T GGG|T|GG|G ||TG|||G|T 

GGIITIIIIT GGTIITITGl G||GG|GGGT I||GG||||G IIGTTGTIIT |GT|||TG|T GGG|T|GG|G I|TG|||G|T 

GGllTllIlT GGTIlTlTGl G||GG|GGGT l||GG||||G ||GTTGT||T |GT|||TG|T GGG|T|GG|G ||TG|||G|T 



gi|21693123|db)|AB076037.1| T||GTGT|G| |G|TTG|TG| T|||GG|TT| |GGGT|||T| |GTTTG|T|T 

AN171 T||GTGT|G| |G|TTG|TGG TT||GG|TT| |GGGT|||T| IGTTTGITIT 

AN169 IllGTGTlGl |G|TTG|TGG TT||GG|TT| |GGGT|||T| |GTTTG|TiT 

AN1S4 lllGfiGllGl |G|TTG|TG| T|||GG|TT| |GGGT|||T| |GTTTG|T|T 



iGGilTG TlllGGTlTG 

iGGiiTG TlllGGTlTG 

iGGilTG TlllGGTlTG 

iGGllTG TlllGGTlTG 




GTTTGlllGl 594 
GTTTGlllGl 600 
GTTTGlllGl 600 
GTTTGlllGl 595 

700 

TllllTGiGG 694 

T||||TG|GG 700 

iGG 700 

tIIIITgIgG 695 



gi|21693123|dbj|AB076037.1| i|||GG|||| 

AN171 T|||GG|||| 

AN169 TlllGGllll 

AN154 TlllGGllll 

gi|21693123|db)|AB076037.1| |GTGG|G|TG 

AN171 |GTGG|G|TG 

AN169 |GTGG|G|TG 

AN154 |GTGG|G|TG 

gl|2ie93123|dbj|AB076037.1| G||G||G|GG 

AN171 GllGllGiGG 

AN169 GllGllGiGG 

AN154 G||G||G|GG 



gi|21 693123|dbi|AB076037. 1 1 
AN171 
AN 169 
AN154 G| 



gl|21693123|dbj|AB076037.1| GG||TGG||T 

AN171 GGllTGGllT 

AN169 GGllTGGllT 

AN154 GG||tGG||T 



G||G|TGG|T TllTlTGTll 
G||G|TGG|T TllTlTGTll 
GllGlTGGiT TllTlTGTll 
GllGlTGGiT TllTlTGTll 



TTG|T|GGTT | 
TTGlTlGGTT | 
TTGlTlGGTT | 
TTGlTlGGTT | 




TiTllGTllT GGlTGlTiTT 

T|T||GT||T GG|TG|T|TT 

TiTllGTllT GGITGItITT 

TiTllGTllT GG|TG|T|TT 

840 

GiTllGllll GIIIIIGTGII 

GITIIGIIII GlIITGTGll 

GiTllGllll GlIITGTGll 

GiTllGllll GlIITGTGll 

I 

IIIEGllllEGil iGlllGllGT 

IIIGIIIIGI iGlllGllGT 

IIHGIIIIGI IGlllGllGT 

llfGlllTGT fGlllGllGT 



iTlllGl llGlTllTGT |TGGGG|TG| |G|T||GG|G |TT|||||GT 794 

IIg|T||TGT |TGGGG|TG| |G|T||GG|G ITTIIIIIGT 800 

G| I|g|t||tgt ItggggItgI |g|t||gg|g IttIIIIIgt soo 

1 IIgItIItgt ItggggItgI |g|t|||g|| |TT|||||g| 795 

_ i|TGTG|T| 894 
iTGGimiT TTflllllll IGITGTGITI 900 

Z ~"~ TTTlllllll iGlTGTGlTl 900 

IGITT IITIIGIIIT TTTlllllll |G|TGTG|T| 895 

960 980 1 000 

GTlTGGlGGi T|GG|GG|GT GG|||G|||T GGGGG||G|G GTG|||G|G| 994 

GTlTGGlGGi T|GG|GG|GT GG|||G|||T GGGGG||G|G GTG|||G|G| 1000 

GT|TGG|GG| T|GG|GG|GT GG|||G|||T GGGGG||G|G GTG|||G|G| 1000 

GTlTGGlGGi T|GG|GG|GT GG|||G|||T GGGGG||G|G GTG|||G|G| 995 




GTTTGGGlTl G|GG|T|||| T|GG||T|TT 
GTTTGGGlTl G|GG|T|||| T|G|||T|TT 
GTTTGGGlTl G|GG|T|||| T|G|||T|TT 
GTTTGGGlTl G|GG|T|||| T|G|||T|TT 



|G iTllGGTTll GG|T|TGG|G G|T|TG|GTG 

|G ITllGGTTll GGlTlTGGiG G|T|TG|GTG GG||G| 

|G ITllGGTTll GG|T|TGG|G G|T|TG|GTG GG||G[ 

|G ITllGGTTll GG|T|TGG|G G|T|TG|GTG GG||G| 



1094 
1100 
T 1100 
1095 



GT|GG|G|TG |G| 
GT|GG|G|TG 
GT|GG|G|TG 
GT|GG|G|TG 




1 200 



IGIOIIIIII lllGllllll IBGTTI 
IGIGIIIIII lllGllllll IIGTTI 
IGIGIIIIII lllGmm IITGTTl 
IGIGTITTTI TTTGllllll TTGTTi 



|G| TlllllGGiT GGTGTlllll TTGGGTTG|T T|||GGG|TG 1194 

|G| TlllllGGiT GGTGTlllll TTGGGTTGIT T|||GGG|TG 1200 

|G| TlllllGGiT GGTGTlllll TTGGGTTG|T T|||GGG|TG 1200 

|G| TlllllGGiT GGTGTlllll TTGGGTTG|T T|||GGG|TG 1195 



gi|21693123|dbj|AB076037.1| { 

AN171 ' 

AN169 I 

AN154 I 



IIGTTTIIII lllGGllGTT TT|G|T|||T ITGTTGTGGI |||T||G||G ||G||T|TTT G|||GGT||| |TG|T|||G| TiGlTlTGiT 1294 

llGBngCIIII lllGGllGTT TTlGllllll |IG|II!GXGG| IXKDGIIgIIG llGTlTlTTT G|||GGT||| ITGItIHgI T|G|T|TG|T 1300 

llGlliim lllGGllGTT TTlGllllll ||6i|G|6G| IBIil|G||G ||GT|T|TTT G|||GGT||| ITGItIHgI T|G|T|TG|T 1300 

■iGlllim lllGGllGTT TflGllllll ||6ilGlGG| |lil||G||G ||GT|T|TTT G|||GGT|T| |TG|T|||G| T|G|T|TG|T 1295 



gi|21693123|dbj|AB076037.1| ||G||GTG|G GG|G|GGG|| GG|G|[I;|[I||!G iTIII |TGG|| 1334 

AN171 ||G||GTG|G GG|G|GGG|| GG|G|Rj|[i!EG IlllllEljlGG GTlTTTTTTi |TGTTG|||| GGTTGGTGGT TGG||G|T|| T|||T|T|TT G|TG|TGG|| 1400 

AN169 ||G||GTG|G GG|G|GGG|| GG|G|||||G i||||li|GG GI|IiElXiI| |TGTTG|||| GGTTGGTGGT TGG||G|T|| T|||T|T|TT G|TG|TGG|| 1400 

AN154 ||G||GTG|G GG|G|GGG|| GG|G|||||G IIIIIIIIGG 6||||||||| liGTTGllll GGTTGGTGGT TGG||G|Ii|| T|||T|T|TT G|TG||G|GG 1395 



gi|21693123|dbj|AB076037.1| 
AN171 
AN 169 
AN154 

gi|21693123|dbi|AB076037.1| 
AN171 
AN169 
AN154 

gi|21 e93123Idbj|AB076037. 1 1 
AN171 
AN 169 
AN154 



IIIGIIG II 

IIIOilG IIOllllGll 

flIlGilG llGilllOll 

llillll IIGIIIIGII 




IITIGIIITG G 

l|6|6||||| |||G|G|i|| illlGGllTG ||T|G|||TG GG|G|T|TTG 
llGlGlllll IliGlGllll IlllGGllTG ||T|G|||TG GG|G|T|TTG 
llGnGlllll linGlGllll IlllGGllTG ||||G|GG|| GGlGlTlTTG 




llTTGGlTTi 1372 
llTTGGlTTi 1500 
llTTGGlTTi 1500 
llTTGGlTTi 1495 



IIGGGIIIIT T|G|-TGG| 

|G|GGG|||| TGT|G|TGG| 

iGlGGGllll TGTlGlTGGi 

IIGGGIIIIT T|G|T|GTG| 



TiGllllfll GGITITITIG GTG||TTTT| G|T||TGG|| TTGTT|||T| ||GTT||G|| G|||T||||T 1470 

TiGllllTll GG|T|TGT|G GTG||TTTT| G|T||T|G|| TTGTT|||T| ||GTT||G|| G|||T||||T 1600 

TiGllllTll GGlTlTGTiG GTG||TTTT| G|T||T|G|| TTGTT|||T| ||GTT||G|| G|||T||||T 1600 

GG|T|TGT|G GTGIITTTTI G|T||TGG|| TTGTTlllTl ||GTT||G|| G|||T||||T 1595 



_ _ IIIIIIIITG GTl 
■G| lllllllllrG GTl 
llGl IIIIIIIBTG GTl 
~ ~ lllllllTTG GTl 



gi|21693123|dbj|AB076037.1| |G|ilTG|T|| HXlXlllTl 

AN171 IGIIIGIIII llllllllll 

AN169 IGIIIGIIII lilililltl 

AN1S4 IGIIIGIIII niililllll 




IGIIIGIGGI Tlllllllll HIIGIIGII 1570 
iGlllGlGGi Tlllllllll IIIIOllGll 1700 
IGIIIGIGGI Tlllllllll illlGllGll 1700 
IGIIIGIGGI TlliniiiGl IlllGllGll 1695 



G| TTlGGTTGTT G||T|G|T|| |T|TG|TT|| | 
G| TTlGGTTGTT G||i|G|l|| |II|IG||]|I 



TiTTGlTlll GG|T|G|G|| TG|TG|||GT 1670 
lIlTGlTlll GGlTlGlGll TG|TG|||GT 1800 



TTlGGTTGTT G||||6|||| |II||G|||| ||||||ll|| l||llG|ni|| GG||||G|G|| TG|TG|||GT 1 
TTlGGTTGTT G||||Gil|G |||||G|||| iGi||||||| |||BG|n||| GG|niG|G|| TG|TG|||GT 1795 



gi|21693123|db)|AB076037.1| | 

AN171 ' 

AN169 I 

AN154 I 



GlHlGlllGG TGGTGGTTGG ||B|G|KT|I GTTGIIgIIg TIIgGGTIIg IIgIgtIIII lllGllllll ||||G|6|6G G|||G|||GG 1770 

III GIIIIGIIIGG TGGTGGTTGG |||||IBT|| GTTG||G||G T||GGGT||G ||G|GT|||| |I|G|||||| IIIIGIGIGG G|||G|||GG 1900 

III GlIllGltGG TGGTGGTTGG ||B||ilFT|| GTTG||G||G T||GGGT||G ||G|GT|||| |||G|||||| ||||G|G|GG GlfllGlllGG 1900 

GlIIITGlTGG TGGTGGTTGG ||T|GTTT|| GTTG||G||G T|TGGGT||G ||G|GT|||| |T|G|T|TT| I|||G|G|GG G|TTGT|TGG 1895 



gi|21693123|dbj|AB076037.1| |G|G|||GT| I|T||GG|TG G||GGG|GG| ||||XGnXiE| ||G||I!E|G|!E 

AN171 IGIIIITGTI llTllGGlTG G||GGG|GG| 

AN169 iGllllTGTl ||T||GG|TG G||GGG|GG| 

AN154 IGllllTGTl I|T||GG|TG G||GGG|GG| 

2 020 

gi|21693123|db)|AB076037.1| HIIIIIgIH G||G|||||G 
AN171 OBIIIIlGlll GllGlllllG 
AN169 illlllGlll GllGlllllG 
AN154 IlllllGlll GinGlllllG 



ITIIIIGIIT GTGlllllTl ||TG|T|TT| |G||G||||T 1870 

IIIIIIGIIT GTGlllllTl IITGITITTI |G|TG|||||! 2000 

GTGlllllTl IITGITITTI |G|TG||||n 2000 

GIGlllllll lllGllliEll iGllGimi 1995 



GHIlGGGlGl lllGlll 1917 

GlillGGGlGl lllGlll 2047 

GKIlGGGlGl lllGlll 2047 

GKIlGGGlGl lllGlll 2042 



Figure 1 Comparative analysis of the zearalenone lactonohydrolase gene sequence in the Trichoderma and Clonostachys isolates 
compared to the complete sequence of the model gene C. rosea AB076037. AN 171, AN 169, AN 154 isolates with identified sequences 
homologous to the zearalenone lactonohydrolase gene, origin - the sequence of the model gene - AB076037. 



Popiel et al. BMC Microbiology 2014, 14:82 
http://www.bionnedcentral.conn/1471 -21 80/1 4/82 



Page 4 of 1 2 



100- 



„ 75- 



c 
o 



50- 



CD 
> 

DC 



25- 




Oh 



1 1 — 

0h5m 2h 



Strain 

— AN154 

- ■ AN 169 

p AN171 



4h 



6h 



Id 



2d 3d 

Time point 



4d 



5d 



6d 



7d 



8d 



9d 



Figure 2 Kinetic reduction of zearalenone during incubation experiments with isolates AN 154 (C rosea), AN 169 (C. catenulatum) and 
AN 171 (7. aggressivum). Experiments were carried out at 25°C, in liquid Czapel<-Dox medium supplemented with yeast extract 
and zearalenone. 



biotransformation ability process following fluctuations in 
transcript levels - see Figure 3). 

Phylogenetic analysis of gene sequences In multiple 
species 

In unsupervised clustering of a/b-hydrolases from mul- 
tiple genomes, newly sequenced homologs of zearale- 
none lactonohydrolase were found to form a stable 
cluster of sequences with previously published homo- 
logues from B. ochroleuca, as well as 2 additional pro- 
teins from M. hrunnea and A. montagnei. 

While phylogenetic reconstruction by maximum likeli- 
hood indicated strong support for a monophyletic clade 
formed by the cluster members (Figure 4), positioning of 
the resulting clade within a/b-hydrolase phylogeny was 
poorly supported and thus remains uncertain. 

Homology modelling and comparative structure analysis 

The created homology models uncovered similarities in 
the active site pocket, as detected by fpocket [15]. In all 
of the modelled structures, the active site pocket is 
strongly hydrophobic under normal conditions - likely 
the catalysis is enabled by allowing access to the active 
site (conformational changes involving cap domain) 
which allows the reaction to proceed by standard 



mechanism involving forming a transient oxyanion hole 
and subsequent cleavage of the lactone ring (Figure 5). 
While homology-based models are likely insufficient for 
elucidation of full sequence of events during substrate 
binding and catalysis (both the variable cap domain e.g. 
[16,17] and surrounding loops [18] are involved in con- 
trolling and fine-tuning substrate access), we were never- 
theless able to ascertain the key functional residues 
involved. 

Our identification of the catalytic triad conflicts with 
the initial proposition of Takahashi-Ando [11] that active 
site is formed by S102-H242-D223 (numeration by align- 
ment in Figure 6). Typically, the nucleophilic attack of 
hydrolase enzyme is facilitated by interaction of histidine 
with acidic residue (third member of catalytic triad). 
This role, according to aU our homology-based models 
cannot be fulfilled by D223 (residue located distantly to 
active site - Figure 7). 

A typically, the third member of catalytic triad appears 
to be El 26 residue, where the side chain is capable of 
interacting with distal nitrogen of catalytic histidine, 
provided conformational changes aUow rotation of the 
glutamate side chain towards histidine (see Figure 5 for 
conformations in modeUed structures). This residue 
is sequentially equivalent (see Figure 7) to catalytic 
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Figure 3 Relative normalized expression (N-fold) of zearalenone lactonohydrolase transcripts during incubation experiments with 
isolates AN 154 (C rosea), AN 169 (C catenulatum) and AN 171 (7. aggressivum). Experiments were carried out at 25°C, in liquid 
Czapek-Dox medium supplemented with yeast extract and zearalenone. 



glutamate residues demonstrated in human epoxide 
hydrolase (PDB:2Y6U, El 53) and epoxide hydrolase from 
Pseudomonas aeruginosa (PDB:3KDA, E169). 

Another possibility is residue D31 - however it ap- 
pears to be nonconserved in Marssonina sequence 
(alanine substitution). Sequencing error cannot be 



completely ruled out in this case, as a single nucleo- 
tide change is sufficient for aspartate to alanine 
substitution in this context. Notably, D31 residue pos- 
ition in relation to the active site histidine favorises 
interactions with proximal imidazole nitrogen (mean 
distance of ca. 2.5 across models) - suggesting 
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Figure 4 Maximum likelihood phylogenetic tree of zearalenone lactonohydrolase homologs from divergent filamentous fungi. 

Bootstrap support is indicated below bifurcations (1000 bootstrap iterations). Tree was based on 245 distinct patterns within a trimmed alignment 
of full length protein sequences (see: Methods section). 
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Figure 5 Superposed structures of template 2XUA (3-oxoadipate 
lactonase; catalytic domain colored in green, cap domain colored 
in yellow) and homology models for zearalenone lactonohydrolase 
homologs from multiple species (see corresponding alignment on 
Figure 6). Coloring is based on RMSD between superposed Ca atoms 
(blue - best, red - worst; gray parts not included in superposition). 

v J 

possible conformational change (freeing the imidazole 
ring) during substrate binding. 

Discussion 

Zearalenone is one of the most dangerous mycotoxins 
produced by fungi belonging to the Fusarium genus. 
Those species are usually severe pathogens of cereals 
and legumes, and may cause Fusarium head blight and 
Fusarium ear rot of corn. These toxins are contributing 
to significant economic losses in livestock production 
causing the disease known as estrogenic syndrome, 
which results in a sterility. Since 1988 [10] it is known 
that among the fungi of Hypocreales order, the mycopar- 
asitic fungus C. rosea have the ability for zearalenone de- 
composition but so far no such properties has been 
described in any species of the Trichoderma genus. 

Selected mycoparasitic Trichoderma and Clonostachys 
isolates were found to be able to reduce significantly 
both the production of zearalenone on medium Czapek- 
Dox broth with Yeast Extract [19] and to detoxify zeara- 
lenone. The three isolates (AN 154, AN 171 - especially 
AN 169) were clearly demonstrated as possible agents 
with verified biotransformation ability {in vitro). This 
finding includes the first demonstration of zearalenone 
lactonohydrolase activity present in a member of Tricho- 
derma genus (AN 171 - T, aggressivum). 



Both gene expression and the ability of isolate AN 171 
(r. aggressivum) to reduce zearalenone levels were con- 
firmed in vitro experiments. The HPLC experiments 
verified that all three isolates are able to decrease toxin 
concentration in a quick and effective manner [11]. The 
RT-PCR analyses further indicate that the expression of 
the zearalenone lactonohydrolase gene is subject to dif- 
ferent modes of regulation in examined isolates. In par- 
ticular, for the isolate AN 171, two hours after the toxin 
administration, a significant increase in the zearalenone 
lactonohydrolase expression is noted, suggesting that in 
T aggressivum the presence of zearalenone in the 
medium directly activates expression of the gene. 

Further study of sequence variation in lactonohydro- 
lase genes is planned, with redesign of PGR markers 
based on sequenced regions and extension into non- 
coding regions of transcript (5'-UTR) [11] using RACE- 
PCR. Subsequent research will also encompass separ- 
ation and identification of end products for detoxifica- 
tion process, as well as isolation of enzyme protein using 
Western blot. 

Previous works have confirmed the existence and 
function of zearalenone - specific lactonase in Clonosta- 
chys sp. (old name of Gliocladium sp.) [9]. The enzyme 
is one of the reasons Clonostachys growth is not inhib- 
ited by zearalenone. We posit that presence of function- 
ing homologues within Trichoderma can also contribute 
to their effective antagonistic activity [19], against 
zearalenone-producing F, culmorum and F, grami- 
nearum (and possibly other resorcyclic acid lactone 
producers). 

Mechanistic features of catalytic site involved in zearale- 
none biotransformation ability are shown to be evolutionar- 
ily old, likely predating the split between Leotiomycetes and 
Sordariomycetes (barring horizontal transfer between fungal 
hosts). While it is unlikely that the exact function of distant 
homologs is the same, the affinity towards large hydropho- 
bic epoxides and conservation of catalytic mechanism (as 
evidenced by active site superposition - Figure 7) are likely. 
Presence of several conserved arginines within the cap do- 
main raises possibility of their involvement in substrate 
binding or orientation (coupled with conformational 
change), analogous to the mechanisms observed previously 
in dienelactone hydrolase [20] and 3-oxoadipate enol lacto- 
nase [16]. Elucidation of the full substrate orientation/ca- 
talysis scenario (including involvement of the glutamate 
and aspartate residues and their spatial conformations dur- 
ing the process) is planned through application of molecu- 
lar dynamics experiments for modelling of the ligand 
binding process. 

Notably, according to previously published work on B. 
ochroleuca enzyme [11] ZEN was rapidly replaced with 
conversion product. The mass of the molecular ion (M + 1) 
corresponding to this product was 293. In our analysis, we 
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Figure 6 Multiple alignment of protein sequences corresponding to: template structure 2XUA (3-oxoadipate lactonase), template 
structure 2Y6U (peroxisomal epoxide hydrolase Lpxl) and lactonase homologs from examined isolates (AN 154, AN 169, AN 171), as 
well as reference sequences from Bionectria ochroleuca (GBK:AB076037), Apiospora montagnei (JGI:58672) and Marsonnina brunnea 
(IVIBIV1_00923 = GBK:EKD21810). Proposed catalytic residues marked with (E128/D31, S102, H242). Residue D223 [11] marked with '!'. 
Secondary structure annotated based on PDB records (2XUA, 2Y6U) and RAPTORX 3-state SSE predictions (a-helix - red, b-sheet - blue). Predicted 
cap domain enclosed in yellow square. 



did not register the corresponding peak, either due to dif- 
ferences in protocol or because of another mechanism of 
zearalenone decomposition. An alternative mechanism of 
zearalenone biotransformation ability was recently demon- 
strated in basidiomycetous yeast Trichosporon mycotoxini- 
vorans [11], which requires a modification of hypothetical 



ZOM-1 intermediate by an unspecified a/b-hydrolase. One 
possibility, testing of which is beyond the scope of current 
work, is that the one-step lactonohydrolase evolved as a 
neofunctionalisation (present within filamentous fungi of 
LeotiomyceteslSordariomycetes orders) of the two-step de- 
toxification mechanism retained by T, mycotoxinivorans. If 
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so, the original mechanism can still exist in select extant 
lineages (within filamentous Ascomycota) in varying degrees 
(dependent on selection pressure towards one-step 
detoxification). 

Conclusions 

Our research shows the first finding of a functional zear- 
alenone lactonohydrolase in mycoparasitic Trichoderma 
aggressivum (an activity earlier characterised in the Clo- 
nostachys rosea strains). Based on the combined screen- 
ing of over ninety isolates of Trichoderma/ Clonostachys 
and in silico investigation of origins of the enzyme activ- 
ity (through phylogeny reconstruction and homology 
modelling) we were able to provide supporting evidence 
for its evolutionary origins, as well as monophyly of 
functional lactonohydrolase homologs in both genera. 
The supporting evidence for presence and activity of 
functional enzyme homologs is based on chemical ana- 
lyses, gene expression patterns, homology models show- 
ing conservation of key structural features and a marked 
reduction of zearalenone content in cultured samples 
(containing both medium and mycelium). 

Methods 

Fungal isolates 

Fungal isolates originated from culture collections of the 
Institute of Plant Genetics (Polish Academy of Sciences, 
Poznan, Poland); Institute of Science of Food Production 
(Bari, Italy; ITEM), Institute of Food Technology 
(Poznan University of Life Sciences, Poznan, Poland), 
Department of Forest Pathology (Poznan University of 
Life Sciences, Poznan, Poland), Research Institute of 
Vegetable Crops, (Skierniewice, Poland) and Rothamsted 
International UK. The isolates were derived from soil. 



compost, wood, cultivated mushroom and cereal grain 
samples. All 98 isolates were identified using both mor- 
phological [21] and molecular methods (ITS 4-5 and tefl 
markers) (Additional file 1: Table SI). 

Isolation of pure cultures 

Fungal isolates investigated in this study were collected 
from pieces of decaying wood, cultivated mushroom 
compost, samples of soil and cereal grain. The samples 
were plated on salt water nutrient agar (SNA) [22] and 
incubated at 20°C for 6 days. Putative Trichoderma and 
Clonostachys colonies were purified on potato dextrose 
agar (PDA, Oxoid). Pure culture were transferred to the 
tubes containing SNA medium and stored at -20°C for 
further study. 

Isolation of DNA 

Mycelium used for DNA extraction was obtained by in- 
oculating Czapek-Dox broth ( Sigma- Aldrich) with Yeast 
Extract (Oxoid) and streptomycin sulphate (50 mg/L"^, 
AppliChem) and after incubation at 25°C for 21 days on 
a rotary shaker (120 rpm). Mycelium was collected on 
filter paper in a Biichner funnel, was held with sterile 
water, frozen at -20°C, and freeze - dried. Total DNA 
was extracted using the CTAB method [23]. The quality 
of DNA was estimated by NanoDrop 2000 UV-vis Spec- 
trophotometer (Thermo Scientific, Wilmington, USA) 
and via Experion Automated Electrophoresis System 
(Bio-Rad, Hercules, CA). 

Primer design 

In the case of C. rosea zearalenone lactonohydrolase, 
previous experiments performed by [9] suggested the use 
of degenerate starters for identification of homologous 



Popiel et al. BMC Microbiology 2014, 14:82 
http://www.bionnedcentral.conn/1471 -21 80/1 4/82 



Page 9 of 1 2 



sequences. In our approach to direct sequencing of ampli- 
fied fragments, degenerate primers gave only non-specific 
products. Because of that seven pairs of primers were de- 
signed on basis of available GenBank homologs (Table 1). 
The primers targeted evenly spread sites along the coding 
sequence (ca. 300 bp estimated product length; estimated 
melting temperature: 60°C). Primer pairs were designed in 
Primer 3 [24] and manually adjusted based on evaluation of 
melting parameters in OligoCalc [25]. 

PGR amplification and sequencing 

The PGR reaction was carried out in a 25 [A reaction mix- 
ture containing the following: 1 [A 50 ng/[A of DNA, 2.5 [A 
10 X PGR buffer (50 mM KGl, 1.5 mM MgGl2, 10 mM 
Tris- HGl, pH8.8, 0.1% TritonX-100), 1.5 \i 100 mM dNTP 
(GH Healthcare), 0.2 \A 100 mM of each primer, 19.35 \A 
MQ H20, 0.25 [A (2 U/[A) DyNAzyme TM II DNA Poly- 
merase (Finnzymes). Amplifications were performed in 
GIOOO thermocycler (BIO RAD, USA) under the following 
conditions: initial denaturation 5 min at 94°G, 35 cycles of 
45 s at 94°G, 45 s at 56°G for all 7 pare primers, 1 min at 
72°G, with the final extension of 10 min at 72°G. Amplifica- 
tion products were separated on 1.5% agarose gel (Invitro- 
gen) in 1 X TBE buffer (0.178 M Tris-borate, 0.178 M boric 
acid, 0.004 M EDTA) and stained with ethidium bromide. 
The 10- (il PGR products were combined with 2 [A of load- 
ing buffer (0.25% bromophenolblue, 30% glycerol). A 100- 
bp DNA LadderPlus (Fermentas) was used as a size stand- 
ard. PGR products were electrophoresed at 3 V cm"^ for 
about 2 h, visualized under UV light and photographed 
(Syngene UV visualizer). The 3 [A PGR products were puri- 
fied with exonuclease I and shrimp alkaline phosphatase ac- 
cording to [26]. Sequencing reactions were prepare using 



Table 1 Sequences of the primers used for amplification 
and sequencing 



Primer name 


Sequences (5-3') 


Estimated product 
length 


LacDP26F 


GAGCCAAGAGAGACCCACAG 




LacDP347R 


TOTGTCCGAATGTCGTOA 


321 


LacDP326F 


GTOAACGACATOGGACAT 




LacDP712R 


AACGTAGTGACCCTGAAGCC 


386 


LacDP693F 


GGOTCAGGGTCACTACG^ 




LacDP903R 


GTATCCTGTCGGGGTAACCG 


210 


LacDP886F 


GTOCCCCGACAGGATACGC 




LacDP1208R 


GAAAGACTCGGTOGTGTCG 


322 


LacDPllBBF 


GCGACACCAACCGAGTCm 




LacDPUOOR 


TACAATATCGCCTGCCCTCT 


212 


LacDP1380F 


GAGAGGGCAGGCGATATOT 




LacDP1695R 


GGGAGCGAGTCAACAACCTA 


315 


LacDP1661F 


AATCTCCGCCATGCTOGG 




LacDP1990R 


GGCTGGTCTCCCGTACAAT 


329 



the ABI Prism BigDye Terminator Cycle Sequencing 
ReadyReaction Kit in 5 [A volume (Applied Biosystems, 
Switzerland). DNA sequencing was performed on an ABI 
PRISM3100 GeneticAnalyzer (USA). 

Sequences were edited and assembled using Chromas 
v.1.43 (Technelysium Pty Ltd). CLUSTAL W [27] and 
MUSCLE [28] were used to align the nucleotide se- 
quences for comparison; the resulting alignments were 
inspected, merged and refined manually. 

RNA isolation and gene expression data analysis 

Mycelium was collected from the Czapek-Dox medium. 
Each sample was weighted on laboratory scales (Sartor- 
ius). Total RNA was purified using RNeasy Plant Mini 
Kit (Qiagen, Hilden, Germany) according to the manu- 
facturers' protocol with the additional DNase digestion 
step. The quality of total RNA was estimated by Nano- 
drop (Thermo Scientific, Wilmington, DE) and via Bioa- 
nalyzer (Bio-Rad, Hercules CA). 

The primer pairs specific to target gene were designed 
using zearalenone lactonohydrolase gene sequences ob- 
tained from T. aggressivum, C, rosea, C, catenulatum 
isolates (Table 2). Analogously to the DNA sequencing 
primers, these were designed with use of Primer 3 [24] 
and their properties were tested using OligoCalc [25] . 

The RT-PCR reaction was conducted using SYBR® 
Green Quantitative RT-qPCR Kit (Sigma- Aldrich). The 
total reaction volume was 25 \A: 12.5 \A SYBR Green 
Taq Ready Mix, 1 \A RNA (< 35 ng), 0.5 \A each primer 
(10 (iM), 0.125 \A reverse transcryptase and 5.125 \A nu- 
clease free water. Gene expression profiles were deter- 
mined through quantitative real-time PCR using a 
CFX96 Touch™ Real-Time PCR Detection System (Bio- 
Rad, Hercules, CA). The reaction was carried using the 
following protocol: initial denaturation 94°C for 2 min, 
followed by 40 cycles at 94°C for 15 s, 61°C for 1 min. 
The melting curve analysis (from 70°C to 95°C) con- 
firmed primer pairs specificity. In the experiment we 
used three biological and two technical replicates to- 
gether with a template-free negative control in each ana- 
lysis of both target and control genes. As a control we 
used mycelium samples cultivated on medium without 
addition of zearalenone. Relative quantification of gene 
expression was done using the 2-AACt method (Bio-Rad, 
Hercules, CA). All data were normalized to |3-tubulin as 
internal control (Real-Time PCR Application Guide, 
Bio-Rad, Hercules CA). 

Mycotoxin chemical analyses 
Sample preparation 

The fungal mycelium was grown in 50 ml Czapek-Dox 
broth (Sigma-Aldrich) with Yeast Extract (Oxoid) for 
9 days at 25°C with rotary shaking at 100 rpm. The zear- 
alenone (Sigma-Aldrich) stock was added after a week of 
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Table 2 The sequences of the primers used for gene 


expression 




Primer 


Sequences (S'-SO 


LACDP723R 


CAAACGTAGTGACCCTGAAGC 


LACDP652F 


CTCGGAGAATGCCAGATG^ 


rtBtubTRICH0R2 


AGCGAATCCGACCATGAAGA 


rtBtubTRICH0F2 


CACCGTCGTOAGCCCTA 



incubation. The initial concentration of ZEA in the li- 
quid cultures was 2 mg/ml The samples (both mycelium 
and medium) were collected before and after addition of 
the toxin. During the first day, the samples were col- 
lected after one minute, two, four and six hours after 
toxin application. In the following days the samples were 
collected once a day at the same time. The collected ma- 
terial was used for both chemical analysis and expression 
profiling. 

Chemicals and reagents 

The zearalenone standard was supplied by Sigma- 
Aldrich-Aldrich (Steinheim, Germany). Acetonitrile and 
methanol (HPLC grade) were purchased from Sigma- 
Aldrich-Aldrich. Potassium chloride was purchased from 
Poch (Gliwice, Poland) and water (HPLC grade) was 
purified with a Millipore system (Billerica, MA, USA). 

Zearalenone analysis 

The samples (lysate containing both medium and myce- 
lia) were filtered through glass microfibre filter (GF/B, 
Whatman). Zearalenone was analysed by the systems 
consisting of: Waters 2695 high-performance liquid 
chromatograph. Waters 2475 Multi X Fluorescence 
Detector and Waters 2996 Photodiode Array Detector. 
Millenium software was used for data processing. The 
excitation wavelength and emission wavelength were set 
to 274 and 440 nm, respectively. The reversed-phase col- 
umn CI 8 (150 mm x 3.9 mm, 4 [im particle. Waters) 
and acetonitrile-water-methanol (46:46:8, v/v/v) as the 
mobile phase at a flow rate 0.5 ml/min were used. Zeara- 
lenone quantification was performed by external calibra- 
tion. The limit of zearalenone detection was 3 (ig/kg. 

The mass spectrometer (Esquire 3000, Bruker Daltonics, 
Bremen, Germany) was operating in the negative ions 
mode with an electrospray ion source (ESI) with the follow- 
ing settings: the source voltage 3860 V, nebulization with 
nitrogen at 30 psi, dry gas flow 9 L min"^, gas temperature 
310°C, skimmer 1: -33 V, MS/MS fragmentation amplitude 
of 1 V ramping within the 40-400% range. Spectra were 
scanned in the mass range of m/z 50-700. The reversed- 
phase column was Alltima CI 8 (150 mm x 2 mm, 3 \im 
particle size) from Alltech. The column was kept at room 
temperature. Three biological and two technical replicates 



were used for each sample. The uninoculated medium with 
added toxin was used as a control. 

Database search and cluster analysis 

The search for zearalenone lactonohydrolase homo- 
logues was conducted on internal, curated MetaSites 
database (Koczyk, unpublished). The dataset consisted 
of combined sequence data from translated GenBank re- 
lease 192 (PLN and BCT divisions) [29], Ensembl/Fungi 
V 16 [30], UniProt/SwissProt [31], PDB [32] and se- 
quences from select, published genomes from JGI/DOE 
MycoCosm [33]. Based on previous BLASTP searches 
for homologs of lactonohydrolase, a single homolog 
from unpublished genome of A. montagnei was included 
in the subsequent analysis. 

The unsupervised cluster analysis was based on the 
subset of proteins detected by 2 iterations of NCBI PSI- 
BLAST [34], on the above-mentioned database clustered 
at 70% protein sequence identity with CD-HIT [35]. 
The zearalenone lactonohydrolase from C. rosea was 
employed as query. The unsupervised clustering of se- 
quences (10728 total) was conducted in CLANS [36], 
using the neural-network based clustering option. 

Multiple alignment and phylogeny reconstruction 

The preliminary alignment of a/b-hydrolases was pre- 
pared with MAFFT [37]. Conserved regions of the align- 
ment were extracted with TrimAl based on 70% gap 
threshold setting [38] . Monophyly of the lactonhydrolase 
cluster within larger context of a/b-hydrolases was then 
assessed with FastTree2 [39] based on LG model (100 
bootstraps) [40]. 

The multiple alignment of zearalenone lactonohydro- 
lase cluster members was prepared using MAFFT-LINSI 
[37], and corrected manually in SeaView [41]. Conserved 
regions of the alignment were extracted with TrimAl 
using 'automatedV setting [38]. Maximum likelihood pa- 
rameters were assessed with ProtTest v3 [42], according 
to Akaike and corrected Akaike information criterions. 
The phylogeny reconstruction for lactonhydrolase ho- 
mologs was conducted in RAxML v 7.3 [43], using 
WAG model of evolution [44], with 1000 bootstrap 
iterations. Template sequence of the oxoadipate enol 
lactonase (PDB:2XUA) was employed as outgroup, in ac- 
cordance with its ESTHER [45] classification in the 
epoxide hydrolase subgroup and its placement in homo- 
logs uncovered by HHpred [46]. 

Visualisation of the phylogenetic tree was prepared 
with ETE2 [47] and custom Python scripts. 

Homology modelling 

Homology modelling was performed with RAPTOR-X 
webserver [48]. Choices of modelling templates were 
checked against HHpred [46] search results for candidate 
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structures in pdb70 (with manual inspection of likely tem- 
plates from epoxide hydrolase superfamily). HHpred was 
accessed via the MPI bioinformatics toolkit portal [49]. 
Visualisation and inspection of all models was conducted 
within PyMol [50]. All structure models are available in 
compressed form in Additional file 2. Multiple alignment 
of zearalenone lactonase homologs is available (in FASTA 
format) in Additional file 3. 

Additional files 



Additional file 1: Table SI. Examined isolates of Trichoderma and 
Clonostachys. 

Additional file 2: Structure models from homology modelling. 
Additional file 3: Multiple alignment of sequences in FASTA format. 
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